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Abstract. NA60 measured dimuon production in p-A and In-In collisions at the 
CERN SPS. This paper presents a high statistics measurement of (j) meson production 
in In-In collisions at 158 AGeV. Both the transverse momentum, rapidity, decay 
angular distributions and the absolute yield were measured as a function of centrality. 
The results are compared to previous measurements in order to shed light on the long 
standing (j) puzzle. In addition, highlights on rj meson production and on the dimuon 
excess below the i /ip mass are presented. 
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1. Introduction 

Lepton pairs production is particularly suitable as a probe of the hot and dense medium 
created in high-energy heavy ion collisions, because they can access several Quark- 
Gluon Plasma (QGP) signatures. Among them is strangeness enhancement [1] that can 
be addressed by measuring the and r] mesons production. The production of thermal 
dileptons [2], that can be detected in the intermediate mass region (IMR) between the 
(j) and the J/tp peaks, is another key observable. 

NA60 is an SPS third generation experiment designed to study muon pairs 
production with high precision and statistics. In the following, after a short description 
of the experimental apparatus, key results on and rj production and on the excess in 
the intermediate mass region will be discussed. 

2. NA60 Detector and data selection 

NA60, a fixed target experiment at the CERN SPS, is mainly dedicated to the detection 
of muon pairs. The apparatus and data selection procedure are fully described in [3]. 
The key detectors are the muon spectrometer (MS) and the vertex tracker (VT). The 
MS detects the muon tracks in the rapidity range 3 < y < 4 , and is preceded by a 
12 A/ thick hadron absorber. The vertex tracker is a Si pixel detector placed between 
the target and the absorber in a 2.5 T dipole field, that tracks the muons, among the 
other tracks, before they suffer multiple scattering and energy loss. The muons detected 
in the MS are matched to the VT tracks both in momentum and angle, thus improving 
the mass resolution and making it possible to measure the distance between the primary 
vertex and the track impact point. 

The sample used for the results presented in this paper was collected with a 
158 A-GeV In beam impinging on a segmented In target. 2.3 ■ 10^ triggers were acquired. 
In order to avoid events with reinteractions of nuclear fragments in the subsequent 
targets, only events with one vertex in the target region were selected. The combinatorial 
background was evaluated by means of the event mixing technique, with an accurancy 
of ~ 1% over the full dimuon mass spectrum. The contribution due to the incorrect 
match between the VT tracks and the muons detected in the MS was estimated either 
with an overlay Monte Carlo simulation or with a mixing technique. The former consists 
in generating a dimuon on top of a real event, while the latter mixes the MS muons 
of a given event with the VT tracks of another event. The two methods agree within 
5%. The final sample consist of 440000 signal muon pairs, with an average signal to 
background ratio of about 1/7. 

3. production in In-In collisions 

NA60 can detect the meson via its decay into both muon and kaon pairs, thus helping 
understanding the so-called puzzle. This takes its origin from the discrepancies 
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observed in the measurements performed by NA50 |1] in the muon channel and by 
NA49 [5] in the kaon channel, both in Pb-Pb collisions at 158 A-GeV. The multiplicity 
measured by NA50 is much higher than the one obtained by NA49. Also the inverse 
slopes Tg// of the transverse mass distributions differ: NA50 finds a value of about 
230 MeV, almost independent on centrality, while NA49 measures values increasing 
with the number of participants Npart from ~250 MeV to ~300 MeV. It has to be 
stressed that the acceptance of NA50 is limited to pt > 1-1 GeV/c, while NA49 is 
dominated by low transverse momenta: in presence of radial flow, the measured Tg// 
will depend on the fit range, being lower at high p^- 

The origin of this discrepancy has been longely discussed. It was suggested that 
kaons may suffer rescattering and absorption in the medium, resulting in a depletion of 
kaon pairs, especially at low pt, that leads to a reduced yield and a hardening of the 
Pt distribution in the hadronic channel [6]. Recent results by CERES [7J both in the 
hadronic and leptonic channels confirm the NA49 results, and thus suggests that there 
is no (f) puzzle. However, the results in dielectrons, being affected by a large error bar, 
are not conclusive. 

NA60 measures the yield in the muon channel by counting the events in the mass 
window 0.98 < M^^ < 1.06 GeV/c^ [8J. The signal/background ratio below the peak 
is ~l/3, integrated over centrality. To account for the continuum under the peak, two 
side windows between 0.88 < M^^ < 0.92 GeV/c^ and 1.12 < M^,^ < 1.16 GeV/c^ are 
subtracted to the window. By performing the measurement in several px, rapidity 
and decay angle intervals, raw differential distributions are extracted. They are then 
corrected for the acceptance using an overlay Monte Carlo tuned iteratively such that the 
resulting distributions reproduce the measured data. The systematic error is evaluated 
by varying the analysis parameters and selections. The measurement is performed for 
5 centrality bins, corresponding to {Npart) = 15, 39, 75, 132, 183. 

The rapidity distribution, integrated over centrality, is fitted with a gaussian, giving 
a width a = 1.13 ± 0.06 ± 0.05. No centrality dependence is observed. The result is in 
agreement with the NA49 results on several colliding systems [5], [9]. The decay angle 
distribution is obtained for several reference systems (Helicity, Collins- S op er, Gottfried- 
Jackson) resulting fiat in all the cases, independent of centrality. 

The transverse mass distributions, shown in fig. [1] (left) for 5 centrality bins, are 
fitted with the function 1/mTdN/dmj- = NQexp{—mT/Teff). The resulting inverse 
slopes are plotted as a function of Npart in fig- [D (right), showing an increase from 
peripheral to mid-central collisions and a saturation for most central collisions. Since in 
presence of radial flow the Te// value may depend on the fit range, in order to compare 
our results with NA49 and NA50, the fit is also performed for pt < 1.6 GeV/c and 
Pt > 1.1 GeV/c. The results are shown in fig. [2] for the two px ranges. At low pt 
(fig. [2], left), the NA60 In-In data show an increase with Npart, in agreement with the 
NA49 results in the same range. At high pt (fig. [21 right) a flattening of the inverse 
slope dependence on centrality is observed, analogously to the NA50 results. However, 
while the difference between the low-pr NA49 results and the high-p^ NA50 results 
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Figure 1. Left panel: <f> transverse momentum distributions in indium-indium 
collisions as a function of centrality. From top to bottom: central to peripheral spectra. 
Right Panel: T slope parameter as a function of the number of participants. 



amounts to ~70 MeV for central collisions, the same difference in the NA60 In-In data 
is of ~15 MeV. While the latter can presumably be due to radial flow, the magnitude 
of the effect in Pb-Pb can hardly be attributed entirely to the same cause. A further 
flattening due to kaon rescattering and absorption may lead to larger Te// values in the 
NA49 data. 



> 350 



300- 



250 



200- 



150- 



> 350; 



300- 



250- 



200- 



150- 



50 100 150 200 250 300 350 400 



50 100 150 200 250 300 350 400 



part 



part 



Figure 2. Centrality dependence of the Te// slope parameter (full circles) compared 
to NA49 (open circles) and to NA50 (open squares). Left panel: fit performed in the 
NA49 range (0 < pt < 1-6 GeV/c). Right panel: fit performed in the NA50 range 
(1.1 <pT < 3 GeV/c). 



The (f) multiplicity was obtained with two different approaches, with largely 
independent systematic errors. In the first, the cross section is directly measured, 
counting the number of produced and the number of incident ions. The average 
multiplicity is then obtained dividing the (f) cross section by the total In-In inelastic 
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Figure 3. {(f)) /Npart in full phase 
space, compared to the NA49 
measurement in several collision 
systems. 



Figure 4. Rcp for the 77 

compared to the one for the cf) 
and the expectation for primary 
production only. 



cross section, estimated with the general formula for two nuclei of mass numbers A 
and B: ainei,AB = vrrg [a^s + _ ^ (^^-1/3 ^ ^-i/sjj ^^^j^ ^j^j^ ^ ^ 3 

(3 = 1.93 which leads to <Jinei,in-in = 4.32 b. Alternatively, the multiplicity can 
be obtained using the J/tp as a calibration process. The ratio between the (p and the 
J/ip multiplicities is given by: (0) / {J/ip) = {N^/A^ x erec,<j,)/{Nj/^p/Aji^ x erec,j/^), 
where N^(^j/^), A^j/^^ and erec.<f>{j/^) are the number of observed (p {J/ip) resonances, the 
acceptance and the reconstruction efficiency respectively. Once corrected for anomalous 
and nuclear absorption, (J/ip) scales with the number of binary collisions Ncou and can 
be written as (J/ip) = {(^nn / n n) ■ Ncoii, where ajj^ and aNN are the J/ip and inelastic 
cross sections in nucleon-nucleon collisions, and Ncou is calculated from Npart using the 
Glauber model. By multiplying this value by the above ratio, one obtains (0). The two 
methods agree within ~ 10%. 

In figl3]the enhancement factor, expressed as the ratio (0) /Npart is plotted as a 
function of Npart, compared to the NA49 results in several collision systems [12], [9]. The 
value in central In-In collisions exceeds the corresponding NA49 one in Pb-Pb. A direct 
comparison with NA50 requires an extrapolation of the NA50 measurements to low pt, 
and is difficult due to the discrepancies in the observed inverse slopes. Two extreme 
cases with T^ff = 220 MeV and T^jf = 300 MeV were considered. The corresponding 
results largely exceed the NA60 values even if the higher T^jf value is considered. 

To conclude, the puzzle is not yet solved. The forthcoming NA60 measurement 
in the KK channel may be decisive to clarify the picture. 

4. ?7 enhancement 

The 7] meson, carrying a strangeness content of ~40%, is expected to be strangeness- 
enhanced in heavy ion collisions. NA60 detects the rj through its Dalitz decay 7] —>■ nnj. 
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In fig. H the Rcp of tlie 77, defined as Rcp{Npart) = iv) /Npart/[{v) /Npart] peripherals 
plotted as a function of Npart, compared to tlie (p. An enliancement of the 17 by a 
factor 1.34 from most peripheral to most central collisions is observed. The values for 
the f] and mesons can be linked through the strangeness under-saturation parameter 
7s. According to [I3], the value of 7^ is ~ 0.7 in Si-Si and ~ 0.9 in Pb-Pb central 
collisions. A value of 0.8 in central In- In collisions is assumed. Since (j) is produced only 
primarily, the dependence on centrality of 7^ can be then extracted from the (p Rcp-, 
according to the scaling: {4>) /N parti/ {4>) /Npart2 = ll{Nparti) hl{Npart2)- The 77 Rcp 
can then be calculated under the hypothesis that 77 is only primary produced: in this 
case Tj/Npart oc QA'~^l{Npart) + 0.6. The result of the calculation is plotted in fig. HJ 
assuming values for 7^ that range from 0.7 to 0.9 in central collisions. It should be noted 
that the results agree with the expectations if only primary production is assumed, in 
contrast to Statistical Models where a large fraction of the rj is produced from secondaries 
without strangeness [H]. For comparison, if the ratio between primary and secondary 
production (the latter coming mainly from the decays of 7ri(1400) and ao which are 
not 7s suppressed) would be 1:1, the expected enhancement from peripheral to central 
collisions would be ~ 1.1. 

5. Results in the intermediate mass region 

The dimuon intermediate mass region (IMR), located between the (j) and the J /ip peaks, 
is expected to be well suited for the search of thermal dimuons, due to its relative 
production yield with respect to the other expected contributions in this region. In 
p-A collisions, the IMR region is described as a superposition of the Drell-Yan and 
open charm sources [15]. The NA38, NA50 and HELIOS-3 experiments observed an 
excess in the IMR with respect to the extrapolation of these sources in S-U and Pb- 
Pb collisions [15], [16]. The mass spectrum of the excess has the same shape as the 
open charm. Possible explanations of this excess are the production of thermal dimuons 
or the enhancement of the charm cross section. In the former case, the excess would 
have a prompt origin, while in the latter the muons coming from charmed meson decays 
would be off- vertex. Due to limited vertex capabilities, it was not possible to distinguish 
between the two scenarios. 

Thanks to its ability to measure the muon offset, NA60 can distinguish, on a 
statistical basis, between prompt and off-vertex dimuons. The resolution of the offset 
distance between the muon tracks and the vertex, estimated using muon pairs from J / %p 
decay, is 37 /xm in the x coordinate (bending plane) and 45 /im in y. Since the offset 
resolution of the matched tracks is affected by multiple scattering and energy loss, a 
weighted offset is used in the analysis: = ^ /\x'^V~^ + Ay'^V^^ + 2/\yI\yV~y' where 
V-^ is the inverse error matrix, accounting for the uncertainties of the vertex fit and 
of the muon kinematics. Ax and Ay are the differences between the vertex coordinates 
and those of the extrapolated muon track z = z^ert- The dimuon offset is defined as 
A^, = ^{Al, + Al,)/2. 
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NA60 observes an excess in the mass region 1.16 < M < 2.56 GeV/c^ with respect 
to the expected sources, that increases from peripheral to central collision faster than 
Npart and Drell-Yan. 

The dimuon offset distribution, shown in figj5] for 1.16 < M < 2.56 GeV/c^, is 
compared to the corresponding distributions for prompt dimuons and open charm, scaled 
with respect to extrapolations from p-A data. It can be noticed that, in order to describe 
the data, the prompt dimuon source must be scaled by a factor 2.43 with respect to 
the expected Drell-Yan, while no enhancement is observed for the open charm, within 
errors. Therefore, the excess can be ascribed to a prompt source. 
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Figure 5. Dimuon 
weiglited offset compared 
to tlie superposition of 
Drell-Yan and open charm 
contributions. 



Figure 6. Inverse slope parameter versus 
dimuon mass. For details, see 1171. 



In order to study the properties of the excess, the expected Drell-Yan and open 
charm contributions were subtracted to the data. The resulting excess mass distribution 
is not reproduced by the Drell-Yan process, being much flatter. The pt spectrum of the 
excess is significantly softer than the Drell-Yan one. Moreover, the excess pt spectra 
extracted from different mass windows in the IMR are clearly different, while for the 
Drell-Yan pairs mass and px spectra factorize, and are thus independent. More details on 
the analysis and results are reported in [3]. The properties of the observed prompt excess, 
being different from the Drell-Yan ones in many observables, suggest an interpretation 
as thermal radiation. 

The inverse trt slope extracted for three mass windows in the IMR, is plotted in 
fig. [6] as a function of the mass, together with the corresponding values obtained in the 
low mass region [T7] . The inverse slopes for the Drell-Yan and open charm contributions 
are shown for comparison. Below 1 GeV/c^, the inverse slope parameters monotonically 
rise with mass from the dimuon threshold where T ~ 180 MeV, up to the nominal 
pole of the p meson, where T ~ 250 MeV. This is followed by a sudden decline to 
~ 190 MeV in the IMR. The initial rise is consistent with the expectation for radial flow 
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of a hadronic source (here tt+tt^ — > p ^ /ip). The jump of ~ 50 MeV observed beyond 
1 GeV/c^, down to a low-flow situation, is extremely hard to reconcile with emission 
sources which continue to be of dominantly hadronic origin. Rather, the sudden loss 
of flow is most naturally explained as a transition to a qualitatively different source, 
implying dominantly early, i.e. partonic processes like qq — > /i/i, for which flow has not 
yet built up. 

6. Summary 

The NA60 experiment measured dimuon production in In-In collisions at the CERN 
SPS. The yield observed in the leptonic channel is higher than the one observed by 
NA49 in the hadronic channel in Pb-Pb collisions. The dependence of the inverse slope 
on pt suggests that the difference between NA49 and NA50 is larger than expected 
from radial flow. The forthcoming NA60 measurement in the — > KK channel may be 
decisive to clarify the picture. 

The r] meson shows an enhancement of a factor 1.34 from peripheral to central 
collisions. This enhancement, compared to the one, is consistent with the expectation 
in the hypothesis that the r] is only primarily produced. 

In the intermediate mass region, an excess of prompt origin is observed. The 
Drell-Yan process can not reproduce the properties of the excess. The value of Te//, 
significantly smaller than in the low mass region, suggests an emission from the partonic 
stage, when the transverse flow has not yet built up. 
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